Excitonic ferromagnetism in the hexaborides 
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A ferromagnet with a small spontaneous moment but with a high Curie temperature can be 
obtained by doping an excitonic insulator made from a spin triplet exciton condensate. Such a con- 
densate can occur in a semimetal with a small overlap or a semiconductor with a small bandgap. We 
propose that it is responsible for the unexpected ferromagnetism in the doped hexaboride material 
Cai-xLa^Be. 



In a recent letter to Nature Young et al^ reported fer- 
romagnetism in Cai_a;Laa;Bg — a material with no par- 
tially filled d- or /-orbitals. Three features are notewor- 
thy, the narrow concentration range of doping (0 ^ x < 
0.01), the small magnetic moment O.OT/is/La), and 
the high Curie temperature (Tc ~ 600 K). The parent 
compound CaBg is a poor conductor which shows a small 
overlap Eq (< 0) between a boron-derived valence band 
and a calcium-derived conduction band at the X-point in 
the Brillouin zonefl The valence and conduction bands 
belong to different irreducible representations (Xg and 
X^, respectively) such that the interband dipole matrix 
element vanishes at the X-point. As a result the dielec- 
tric constant (k) and the effective masses {ruh and mg) 
remain fiiute as Eq 0, which, as pointed by Halpcrin 
and Rice,a favors formation of an excitonic condensate. 

An excitonic instability for a semimetal or a semicon- 
ductor with small values of Eq was first proposed by 
Keldysh and KopaevQ and des Cloizeauxla and later ex- 
amined, in detail by a number of authors in the mid- 
sixties. □ The effective mass Hamiltonian for a gas of elec- 
trons and holes, attracting each other through a Coulomb 
interaction screened by the background dielectric con- 
stant K, is treated in the Hartree-Fock approximation. 
The characteristic energy scale is set by the exciton bind- 
ing energy i?cx = Ai*'«~^Ry, where ^* is the reduced mass 
(/i*^^ = TO^^+m^T^) in electron mass units. In CaBg typ- 
ical values for the effective mass tensors are m\ — 2.17, 

771^ = 0.206 (valence band) and ml = 0.504, = 0.212 
(conduction band) from our band structure calculations 
and At 5 from optical measurements (L. Degiorgi, pri- 
vate communication) giving E^x ~ 0.08 eV. A triplet 
exciton condensate is, favored over a singlet one by the 
Coulomb interaction.cl In the present case of a direct gap 
the broken symmetry will give rise to a local spin po- 
larization in each unit cell but no net magnetization at 
stoichiometry. As first shown by Volkov et al.,U extra 
electrons doped into the system are distributed asym- 
metrically between two spin projections in order to keep 
most favorable pairing conditions for at least one spin 
species of exciton. As a result, the excitonic insulator 
becomes ferromagnetic with the Curie temperature set 
by the energy scale of the primary excitonic order pa- 
rameter. 



The theory of a Bose condensate of loosely bound 
electron-hole pairs has a close formal resemblance to the 
description of Cooper pairs in superconductors. We con- 
sider a valence band maximum and a conduction band 
minimum both situated at k = 0. In the normal phase, 
without excitonic condensate, the electrons states are 
classified by their wavevectors and band numbers with 
creation operators and 6j^^ for valence and conduc- 
tion bands, respectively. The band quantum number 
ceases to exist in the excitonic phase because of the for- 
mation of electron-hole pairs. This change in the ground- 
state is reflected in a new Hartree-Fock average (al^^b^^,), 
which becomes the primary order parameter of the exci- 
tonic phase. Quasiparticles in two new bands are related 
to the old Bloch states via a canonical transformation 



(1) 



where the unitary matrix Afg-o-' describes a spin sym- 
metry of the excitonic condensate. The singlet state of 
excitons corresponds to Ma-a' = ^aa'i whereas a triplet 
exciton condensate has Ma^i = niu]^^, , (t* being the 
Pauli matrices 




FIG. 1. The structure of energy spectrum in a semimetal 
(left) and in an excitonic insulator (right) for one spin po- 
larization. States in shaded volumes are occupied by doped 
electrons. 

We keep only the intraband exchange part of the 
Coulomb potential between electrons and holes, V'(k), 
which has smallest momentura, transfer and is respon- 
sible for electron-hole bindingfl The formation of exci- 
tons opens a gap in the quasiparticle spectrum Ako- = 
^^^(k — q)(akcr^kCT') and the system transforms into 
an excitonic insulator see Fig. 1. Two new bands have 



energy dispersions given by E, 
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with ?7k = ^{e^ + e^) and = 
determined from the equation 



2 V^k 



The gap is 
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where n^^ 



akcr) and n^^ = (/S^^/^kcr) are occupa- 
tion numbers for a- and /3-branches of spectrum. The 
neglected Coulomb terms with a large momentum trans- 
fer select the triplet state with an arbitrary orientation 
of polarization vector nfl 

We start with the simplest weak-coupling variant: 
isotropic bands with masses to/j = rTig = m so that 
eJJ = — — So — ^ — Eq ^ —Eg/2. In the stoichio- 
metric case the number of electrons equals the number of 
holes and /j, = 0. For a significant overlap there is a strong 
screening of the Coulomb interaction by free carriers and 
we can set y(k) « Vq and, hence, neglect a momentum 
dependence of the excitonic gap. At zero temperature the 
lower band is occupied (ng = 1), while the upper band is 
empty (n^ = 0). Then, the spin-independent solution of 
the gap equation is Aq = 2ec exp[— (Vb-/Vo)~^], where Sc is 
a cut-off around the Fermi surface and A^'o = rnkp / {^tt"^). 
The energy gain of the excitonic condensate over the nor- 
mal phase per spin direction is E[Q\ — — iiVoAQ. 

Electron-hole asymmetry created by doping leads to 
partial occupation of the upper /3-band (Fig. 1), which 
according to Eq. (||) gradually destroys the excitonic con- 
densate. We relate the number of doped carries Ud with 
a change of the chemical potential in the normal state: 
Ud = 4:N()fid- Then, the zero temperature excitonic gap 
in the unpolarized state is 



A^ =Ao(Ao-2|/Xd|) . 



(3) 



The excitonic condensate becomes energetically unfavor- 
able for Hd > Ao/2. The ground state energy per spin 
direction (relative to the stoichiometric normal state) is 



(4) 



Polarization is included in this picture by allowing differ- 
ent chemical potentials for spin-up and spin-down quasi- 
particles: ^^(1 + p) and /id(l — p). A quick check 
shows that the energy of a polarized excitonic insula- 
tor E[iidi^ + E[fj,d{l- ~ p)] has a minimum value for 
p = 1, i.e. for complete polarization. All doped electrons 
have parallel spins leaving stoichiometric conditions for 
carriers with opposite spins. This gains pairing energy 
relative to the unpolarized state value 2E[fid\ and com- 
pensates for an increase in kinetic energy. The two gaps 
are nonzero in the following doping ranges 



Ax = ^Ao(Ao - 4/Xd) for < < Ao/4 ; 



At = Ao 



for < < Ao/2 



(5) 



Note, that at intermediate concentrations electrons and 
holes are paired only for one spin direction. At the crit- 
ical point /ij^ = Ao/2 the ground state energy becomes 



degenerate with respect to the degree of polarization p 
and after a jump the system returns back to the normal 
unpolarized state. 

In the case of CaBg the excitonic gap, as estimated 
from Tq is of the same order of magnitude or even 
larger than the overlap \Eg\ — 2eo. Under this cir- 
cumstances the weak-coupling approach does not work 
quantitatively. In addition there is a significant mass 
anisotropy in both bands. We, therefore, have solved 
the gap equation (J2|) numerically using effective masses 
from our band structure calculations and the gap value 
Ao — 0.07 eV as suggested by optical measurements on 
the 5%-doped compound (L. Degiorgi, private commu- 
nication). The band overlap remains a largely unknown 
parameter, which we take to be Eq = —0.2 eV. In ad- 
dition, we again neglect momentum dependences of Ak 
and V^(k), i.e. assume a strong screening with a cut-off 
energy Sc = \Eg\/2. The variation of the two gaps under 
doping is shown in Fig. 2. For small donor concentrations 
the excitonic insulator is again completely polarized. The 
major difference with the isotropic weak-coupling model 
is a new concentration range, where the upper gap grad- 
ually decreases and doped electrons are only partially po- 
larized. Note, that ferromagnetic polarization disappears 
before the transition to the normal metallic state. 
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FIG. 2. The doping variation of the gaps for spin-up and 
spin-down quasiparticles (dashed lines) and of the saturated 
ferromagnetic magnetization (solid line) in Cai_a;Laa;B6 
at T = 0. 

The ferromagnetic moment of the polarized doped ex- 
citonic insulator is = X^kcr '^('^ka + '^kir)' where the 
occupation numbers are to be expressed in terms of the 
eigenstates (^. The exchange electron- hole attraction 
does not fix the relative orientation of the excitonic polar- 
ization vector (n) and the ferromagnetic moment (M^). 
In order to determine the angle between n and M we have 
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to invoke extra Coulomb terms. The principal role in lift- 
ing this degeneracy is played by a direct Hartree contribu- 
tion, which is written in terms of the Fourier harmonics of 
the charge density pc as Vdir = 5 Z^g.g' ^^g.G' {0)p*G,Pc- 
In the case of weak doping, we obtain from (0) the fol- 
lowing expression 

^G = E^-G(k,k)uk.^'k.(M,, + A/;jK,-nfj , (6) 

kcr 

where F^''f.{'k,]i.) are the form-factors defined in Ref. M. 
In order to minimize the direct Coulomb energy charge 
modulations pc must vanish. For an undoped system, 
when occupation numbers and coefficients Ukcr and t;kcr 
are spin independent, this requirement leads to TvM = 0, 
i.e. Vdir favors triplet states over the singlet one. For a 
doped excitonic insulator, n^^ — n^^~ ^ n^^ — n^^ and 

Pg becomes nonzero for M = o"^, whereas it still van- 
ishes for M = (T^ or CT*'. We conclude, therefore, that 
the perpendicular orientation of n and M has the lowest 
energy. 

For transverse polarization of triplet excitons an in- 
duced ferromagnetic moment is 

kcr 

In this case each of the quasiparticles @) consists of 
an electron and a hole with opposite spins and, there- 
fore, carries only a small fraction of the Bohr magneton. 
In addition AI^ has an important contribution from a 
change of the vacuum for spin-up and spin-down quasi- 
particles. The contribution vanishes exactly for the con- 
stant density of states. We used the above expression for 
to calculate the variation of the ferromagnetic mo- 
ment in CaBg under doping, which is presented in Fig. 2. 
The magnetization depends on electron-hole asymmetry 
in the density of states. Smaller overlap leads generally 
to more asymmetry and, consequently, to an increasing 
value of M^. 

For the above choice of parameters we calculated from 
Eq. (||) concentration dependence of the excitonic tran- 
sition temperature for Cai-ajLajjBg and the stability 
boundaries of the ferromagnetic polarization, which are 
shown in Fig. 3. The Curie temperature grows with dop- 
ing and the maximal Tq is achieved above the concentra- 
tion at the peak value for the saturation magnetic mo- 
ment. The direct transition from the normal semimetallic 
state into a polarized excitonic phase, which exists for a 
certain range of doping, will split into two seperate tran- 
sitions when the degeneracy lifting Coulomb terms are 
included in Eq. (||). 

Excitonic ferromagnetism offers a natural explanation 
of weak ferromagnetism of Cai_a;Laa;B6. Recognizing 
that the dimensionless parameter rs (ratio of interparti- 
cle spacing pp to the effective Bohr radius a^) has a mod- 
erate value,EI leads us to conclusion that our explanation 
is more plausible than the original interpretation based 



on a ferromagnetism of a low-density electron gas.Bfl Cal- 
culations within the simplest model of an excitonic insu- 
lator yield values of the critical donor concentration, the 
weak ferromagnetic moment, and the Curie temperature 
close to the experiment. 
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FIG. 3. Phase diagram of the doped excitonic insulator. 
Shaded area shows stability region of the excitonic ferromag- 
net. 

We didn't attempt detailed theoretical predictions 
mainly because of the limited experimental data available 
at present. A major uncertainty is the value of the band 
overlap Eq in CaBg (and in isoelectronic SrBg). Our 
LDA calculations gave a small overlap = —0.2 eV, 
whereas the modified LDA+U methocO gave a small 
band gap Eq — 0.15 eV. In the original theory of the ex- 
citonic insulator Eq was considered to be a band param- 
eter that varied continuously through the value Eq = 0. 
Subsequently, the discovery of the electron-hole liquid 
showecO that a first order transition between states with 
a substantial band gap {Eq > 0) and overlap {Eq < 0) 
should occur and that smaller values of \Ec\ ^ -Eex lie 
in the unphysical intermediate region. At present there 
are insufficient experiments to reach a firm conclusion 
on this point for the hexaboridcs. In addition theoreti- 
cal questions such as the role of the multiple X-points in 
the Brillouin zone and of optical phonons require further 
study. 

Excitonic ferromagnetism offers a qualitative explana- 
tion for the three key features of the ferromagnetism of 
Cai_2;Laa;B6 enumerated in the beginning. Experimen- 
tal predictions which will be a decisive test of this model 
are optical edges at the band gaps 2Af and 2Aj in the 
infrared optical conductivity of doped hexaborides and 
also lowlying magnetic excitations in the stoichiometric 
compounds CaBg and SrBg, which are the spin- waves of 
the triplet exciton condensate. 
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